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Using data from the FOCUS (E831) experiment at Fermilab, we present a 
new measurement of the weak decay-asymmetry parameter in A+ —> 
Att'*' decay. Comparing particle with antiparticle decays, we obtain the 
first measurement of the CP violation parameter A = ——. We obtain 

^ “Ac 

OAc = —0.78 ± 0.16 ± 0.13 and A = —0.07 ± 0.19 ± 0.12 where errors are 
statistical and systematic. 


1 Introduction 

Parity violation in the weak decay of a spin ^ hyperon into a spin ^ baryon 
and a pseudoscalar meson is well known. An example is A —pvr” [1], where 
the angular distribution of the proton in the A rest frame is given by 


dW 



( 1 ) 


dcosO 
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where P is the polarization of A along the z direction, 9 is the polar angle and 
oa is the weak decay-asymmetry parameter. The latter is defined as 

27^e(A;A_) 

| A +|2 + |/ l _|2 

where A+ and are the parity-even and parity-odd decay amplitudes. In 
a non-relativistic picture, they correspond to the L = 0 and L = 1 orbital 
angular momenta of the proton-vr system respectively. 


( 2 ) 


Eqn. (1) can be generalized for the case of a baryon double decay chain like 
Att'*' with A —> pTi~ in which each baryon weak decay is a | | -|-0“ 

process. If the A+ is produced unpolarized, Eqn. (1) for the A decay becomes 

dW 1 

- -- =-(1+ aAc«ACOs0) (3) 

d cos 9 2 

where oa^ is the weak-decay asymmetry parameter of the A+ —> ATr"*" process 
and 9 is the helicity angle of the proton in the A rest frame, i.e. it is the 
supplement of the angle (i.e. vr minus the angle) in the A rest frame between 
the proton and the direction of the A+. For this decay of the A+, this is the 
same as the angle, in the A rest frame, between the proton and the 7r+ decay 
product of the A+, and the supplement of this angle is used in this analysis. 
The distribution of the helicity angle, 9, determines the product oAc^a, and, 
since oa is known, (in this analysis we assume oa = —cnj = 0.642 ±0.013 [2]), 
one can extract oa^- 


If CP were conserved exactly, of the A+ —> A7r+ process would be the 
negative of a— in A“ —Att” decay. Just as with the A decay, there is the 
possibility of CP violation in A+ decay (a weak decay with possible final state 
interactions). In this paper we 1) measure oa^ and a— separately, 2) use them 
for the first measurement of the CP asymmetry parameter 


^ ^ Q^a, ± QaJ 
“Ac - “at 


( 4 ) 


and 3) having established that the difference is negligible within our errors, 
combine the data for particle and antiparticle decays to obtain the best value of 
the asymmetry parameter. In this analysis the A+ is assumed to be unpolarized 
and thus the A longitudinal polarization (polarization along the A momentum 
direction in the A+ rest frame) is ctAc- As discussed later, if the polarization is 
the maximum allowed in our data, the effect on our measurement is negligible. 


FOCUS is a charm photoproduction experiment, an upgraded version of E687 [3], 
which collected data during the 1996-97 fixed target run at Fermilab. Elec¬ 
tron and positron beams obtained from the 800 GeV/c Tevatron proton beam 
produce, by means of bremsstrahlung, a photon beam (with typically 300 GeV 
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endpoint energy) which interacts with a segmented BeO target [4] . The mean 
photon energy for triggered events is ~ 180 GeV. A system of three multicell 
threshold Cerenkov counters is used to perform the charged particle identi- 
hcation, separating kaons from pions up to 60 GeV/c of momentum. Two 
systems of silicon microvertex detectors are used to track particles: the hrst 
system consists of 4 planes of microstrips interleaved with the experimental 
target [3,5] and the second system consists of 12 planes of microstrips lo¬ 
cated downstream of the target. These detectors provide high resolution in 
the transverse plane (approximately 9 /rm on the track position), allowing 
the identihcation and separation of charm production and decay vertices. The 
charged particle momentum is determined by measuring the deflections in two 
magnets of opposite polarity through hve stations of multiwire proportional 
chambers. 


2 Analysis of the decay mode A+ —> A7r+ -|- charge conjugate. 


Unless explicitly stated otherwise, both the particle and its charge conjugate 
are implied. The selection of Ajf events begins with the identihcation of A 
candidates on the basis of vertexing and loose Cerenkov cuts on the pion and 
proton decay particles as described in detail elsewhere [6] . The tracks of these 
charged daughters are used to form the decay vertex of the A and to determine 
its hight direction and momentum. The momentum of the A is used together 
with the momentum of a tt’*' (or with the momentum of a 7r“ in the case of a 
A) to hnd the A+ (A“ ) vertex and its momentum. We call this pion tt. The 
A’s used in this analysis can be divided into two categories: the ~ 10% which 
decayed before the microstrip detector and the large majority that decayed 
after the hrst microstrip plane. For the former, both the A decay vertex and 
the momentum vector are measured precisely by the FOGUS spectrometer; 
for the latter only the momentum vector is measured precisely. Thus different 
algorithms are used to hnd the A+ decay vertex and the production vertex 
in these interactions. For the A’s decaying upstream of the the microstrips, 
both the vertex position and the direction of the A are used in combination 
with a fully reconstructed tt track. The conhdence level that the tt and the 
A originated from the same Aj" decay vertex is computed and required to be 
> 1%. The A+ candidate decay vertex and momentum are calculated and used 
as the initial track for a candidate driven vertex algorithm [3] . The momentum 
and position of the resultant A+ candidate are used as a seed track to intersect 
the other reconstructed tracks and to search for a production vertex. The 
conhdence level of this vertex is required to be greater than 1%. 

For the A’s decaying downstream of the hrst plane of the microstrip, we use 
the momentum information from the A decay and the silicon track of tt to 
form a candidate A+ momentum vector. This vector and the tt track form a 
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plane in which the momentnm vector of the candidate A+ lies. The transverse 
distance from this plane is calculated for candidate production vertices which 
are formed by at least two other silicon tracks and the conhdence level that the 
vertex lies in the plane is computed. The A+ trajectory is forced to originate 
at this production vertex and the conhdence level that it verticizes with tt is 
calculated. If this conhdence level is greater than 2%, this vertex is taken as 
the decay vertex of the A+ and used with the A+ momentum as the seed track 
for a candidate driven search for the production vertex. The conhdence level 
of this vertex is required to be greater than 1%. 

For both types of A’s, once the production vertex is determined, the distance / 
between the production and decay (A+ decay) vertices and its error ai are com¬ 
puted. The quantity l/ai is an unbiased measure of the signihcance of detach¬ 
ment between the production and decay vertices. This is the most important 
variable for separating charm events from non-charm prompt backgrounds. 
In this analysis a cut of //cx; > 3 has been imposed on the selected events. 
Signal quality is further enhanced by imposing tighter Cerenkov identihcation 
cuts on the proton decay product of the A (to eliminate the contamination 
of Kg decaying into two pions) and tt. The Cerenkov identihcation cuts used 
in FOCUS are based on likelihood ratios between the various particle iden¬ 
tihcation hypotheses. These likelihoods are computed for a given track from 
the observed bring response (on or oh) of all the cells that are within the 
track’s (/3 = 1) Cerenkov cone for each of the three Cerenkov counters. The 
product of all firing probabilities for all the cells within the three Cerenkov 
cones produces a y^-like variable Wi = —2 In(likehhood) where i ranges over 
the electron, pion, kaon and proton hypotheses [7]. The proton track is re¬ 
quired to have IdC — kFp greater than 4 and Wk — Wp greater than 0, whereas 
the TT track is required to be separated by less than 6 units from the best 
hypothesis, that is pion consistency, W m^r, — UC- < 6. Signal quality is fur¬ 
ther enhanced by requiring that the momentum of the A+ be greater than 40 
GeV/c, the position of the production vertex to lie within the target hducial 
volume, the proper lifetime of the candidate A+ to be less than 5 times the 
nominal A+ lifetime, the transverse momentum of the A+ with respect to the 
average beam direction be greater than 0.2 GeV/c and hnally by requiring 
that I cos a; I < 0.8, where u is the angle between the A+ hight direction and 
the TT hight direction evaluated in the A+ rest frame. 

Using the set of selection cuts described, we obtain the invariant mass distri¬ 
bution for Avr” shown in Fig. 1 with 776 ± 55 events in the particle channel 
and 637 ± 34 events in the antiparticle channel. The events in these plots are 
further divided into four slices of the cosine of the A decay angle, and the 
numbers of Avr decays of the A+ above background in each of the four slices 
constitute the histogram of events vs. A decay angle. In order to do the hts 
for the numbers of decays, three additional features must be recognized: 1) 
the change in slope at 2.15 GeV/c^, 2) the A+ —»■ rehection just below 
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the A+ peak, and 3) the smooth background. 


The decay mode A+ —> A7r~^7r^ where the tt^ is not detected contributes to the 
rise below 2.15 GeV/c^. Examination of Monte Carlo generated A+ —> 
events verihed that this reflection becomes negligible above 2.15 GeV/c^ and 
motivates the choice of 2.15 GeV/c^ as the lower limit of the mass region used 
for this analysis. 


The bump from 2.15 to 2.25 GeV/c^ is the decay A+ —> where the 7 

from > Ay is not seen. Including the region provides a check on the 
background parameterization but requires that the angular distributions of 
the mode must be used in the mass and efficiency calculations. Because 
of the electromagnetic decay [8,9] in this chain, the decay is isotropic and, 
integrated over all decay angles, the helicity angle distribution of the proton 
from the A is flat regardless of the initial polarization of the A+ and the value 
of aA^^so-TT (the weak decay-asymmetry parameter of the A+ —> decay). 

However, in this analysis, the acceptance, efficiency, and calculation of the 
cosine of the A decay angle from experimental data may depend upon the E° 
decay angles. Moreover, the helicity angle is miscalculated as the angle, in the 
A rest frame, between the proton and the negative of the direction of the pion tt 
rather than the negative of the direction of the photon. Thus the distribution 
of observed proton helicity angles in this mass region could be distributed 
unevenly among the cos 6 ' slices with the distribution depending on 
These effects have been included in this analysis by generating Monte Garlo 
A+ ^ EO events and reconstructing them as A+ —> A 7 r+ events. The events 
are generated assuming the A+ is produced unpolarized and for a set of values 
of aAc^so-TT which covers the physical range. The differences among the four 
regions of A decay angle are small but have been included in the analysis. 


The background shape is obtained from the wrong sign Avr” (he. the pion from 
the decay vertex and the nucleon from the A have the opposite charge) mass 
distribution. This mass distribution is htted with a second degree polynomial 
and the polynomial is scaled by a multiplicative factor in the subsequent hts. 
The cuts used for this background distribution are the same as for the signal 
except for the additional cut cos 9 < 0.7 to remove events in which the 7 r“ 
(tt'*') track associated by the tracker with a proton to form a A (A) vertex was 
actually a random track, and the true 7 r“ (vr^) track was actually misidentihed 
as TT. No signal is expected or observed in the A+ mass region. For the A+ —> 
Avr'*' signal, two Gaussians with common mean and different resolutions are 
used. These resolutions and relative amplitudes have been obtained with a 
Monte Garlo simulation and normalized so that the sum has unit area for 
each cos 9 slice. In this way the Monte Garlo shape of the signal is preserved, 
and the mass £t varies only an overall multiplicative factor. 
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3 Fitting procedure and extraction of and A 


In order to extract and A , the data sample is divided into particle and 
antiparticle subsamples and each subsample is further divided into four equal 
slices of cos^, spanning the range from —1 to +1. The Avr mass plots of 
the resulting eight subsamples are shown in Figs. 2 and 3. The number of 
A+ —> A7r+ events above background in these plots is the basic data needed to 
measure the slope and thus the product aAOAc^AOTr which is called a in this 
discussion of the £t procedure. These numbers are predicted by Eqn. (5) below 
in a simultaneous fit to all four slices. More specihcally, if Uj is the integral of 
the Gaussians signal term for the i**^ slice, then: 

PCOS Oup^i 

Hi = Nsige{6i, a) / -(1 + Q;cos6')dcos6' (5) 

2coseiow,i 2 

where cos6*up,i and cos6'iow,i are the upper and lower limits of the cos 9 
slice, a = aA^aA is a £t parameter, e{9i,a) is the acceptance-efficiency factor 
calculated with a Monte Carlo for the i*'^ cos 6 slice and is the total number 
of A+ —> Avr^ events produced in the FOCUS experiment (before correcting 
for apparatus acceptance). Since e depends upon a, it has been calculated with 
Monte Carlo for a set of values of a and interpolation provides the continuous 
function needed for the ht. 


The amplitude of the reflection is parameterized as the branching ratio: 


BR5]o 


F(A, ^ SOvr) 
F(A, ^ Avr) 


Arefl 


( 6 ) 


and iVi-efi multiplies the Monte Carlo shape of the reflection. The Monte Carlo 
shape includes the efficiency and is normalized to one event before acceptance 
corrections. Since the acceptance corrections depend slightly on ccAc^sOtt; a 
table is constructed and linear interpolation is used to permit calculations 
with different asymmetry parameters. 


The background is parameterized by a single parameter Ki for each slice which 
multiplies the wrong-sign background. 

In summary, the free fit parameters are : a, oAc^sOtt , Aa^ = J2i=i BR^o and 
the four iFj’s. Finally, the ht is performed maximizing the likelihood function 


4 mass bins 


£=n n e "* 

1=1 j=i 



( 7 ) 


where the sum of the three terms just described is = predicted events in 
the mass bin of the mass plot {i = 1, ...,4, corresponding to each cos6' 
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slice) and sj^i = the number of experimental events falling in the mass bin 
of the i**^ mass plot. 

The £t results are: = —0.73 ± 0.21 (with a x^/d-o.f. of 1.1) and a— = 

0.84 ± 0.22 (with a x^/d-o.f. of 1.0) where the errors are statistical. These 
numbers lead to ^ = —0.07±0.19. Assuming CP invariance in this decay and 
taking a weighted average (after changing the sign of o.—), one obtains = 
—0.79±0.15. The relative branching ratios are BR^o = 1.14±0.15 for particles 
and BRso = 0.81 ± 0.16 for antiparticles, leading to a weighted average of 
BRso = 0.98 ± 0.11 The determined by the £t were at the edge of 

the physical range both for the A+ (aAc^so-Tr = ~1) and for the A“ (a^ 
sample, and with a very large error of 50%. In Figs. 2 and 3 the results of the 
£t are shown for the particle and antiparticle sample respectively, with the 
contribution of the A+ signal and the A+ —> identihed. The amplitudes 

of the signal and of the reflection have been computed from parameters of 
the global ht over all four slices. A further visualization of the quality of hts 
are shown in Figs. 6 and 7. The number of signal events in each slice has 
been computed by subtracting the background and reflection from the data 
histograms in the region between 2.255 and 2.315 GeV/c^. The amplitudes of 
the background and reflection were taken from the global ht. The resulting 
numbers are then corrected for acceptance and normalized for comparison 
with the function 


2/i = 1 ± aA^ttACOs^i (8) 

where yi are the subtracted and corrected data points in the cos 9 bin and 
cos6*j is the value of cos 9 in the middle of the i*'^ bin. The solid line is the 
function in (8) with the value of returned by the ht. 


4 Systematic checks 


Many checks have been performed in order to assess the systematic errors on 
ctA^ and A . 

In this analysis we assume the A+ is produced unpolarized. We think this is a 
safe assumption for the FOCUS experiment. In fact if the A+ were polarized 
the A polarization would be larger than ccAc in one hemisphere about the 
A+ polarization and smaller in the other hemisphere, but the polarization 
averaged over all A+ decay angles remains q;Ac- Even a large, asymmetric error 
in the acceptance correction and a polarization as large as 0.5 would contribute 
negligibly to the systematic error. Preliminary polarization measurements [10] 
with this data indicate a much smaller polarization. Therefore the effect on 
our measurement of a possible A+ polarization is neglible. 


The fitting conditions and the cuts were changed in a reasonable manner, on 
the whole data set. The and A values obtained by these variants are all a 
priori equally likely, therefore this uncertainty can be estimated by the r.m.s. 
of the measurements [11]. 


As explained earlier, the background shape of each Avr^ mass plot was fit 
with a second degree polynomial extracted from the wrong sign sample and 
only an overall scale factor was varied in the mass ht. To check how critical 
this assumption is, a mass ht was performed in which all twelve polynomial 
coefficients (three coefficients times four cos 9 bins) were ht free parameters 
instead of taking the shape from the wrong sign sample. The results of the 
ht performed in this way are A = 0.06 ± 0.17 and = —0.84 ± 0.14 when 
one assumes CP invariance. These results are in agreement with those of the 
standard analysis and show that the result is not sensitive to the way the 
background is ht. The following further systematic checks were performed: 


- using a 2 MeV/ A binning in the mass plots instead of the 5 MeV/ A binning 
of the standard analysis; 

- dividing the cos 9 range in 5 equally wide intervals instead of 4 intervals; 

- varying the most relevant selection cuts in order to check how well the 
Monte Carlo reproduces the experimental apparatus, namely : variation of 
the I/a cut; variation of the cut on the transverse (with respect to the beam 
direction) momentum of the A+; variation of the cut on | cosu;|; sharpening 
of the picon cut; variation of the cut on the proper time of the A+ candidate; 
variation of the Cerenkov cut on the proton. 

- Another check was performed to assess the possible systematic error due 
to errors in the correcting function e{9i,Q) caused both by the interpola¬ 
tion method and by an intrinsic failure of the Monte Carlo to simulate the 
experimental apparatus. The correction functions that differ the most are 
those corresponding to the extreme values of ctAc Q:a, (or a^aj) namely 

®A = —0.9 and oa = 0 (and correspondingly for the A^ sample) 
as shown by Figs. 8 and 9. The ht was performed once imposing cta^ cya = 
a-a-^ = —0.9 (no interpolation) and then imposing oa^ oa = cty^^A — 0- 
These results were taken as an indication of the magnitude of a possible 
systematic error. 


The results for oa^ and A for all these systematic checks are summarized in 
Figs. 10 and 11. As explained previously the r.m.s. of all those results is taken 
as the estimate of the systematic error. For oa^ this is ±0.12 and for A it is 
±0.13 
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Experiment 


FOCUS (this result) 

CLEO II [12] 

ARGUS [13] 

CLEO [14] 

-0.78 ±0.16 ±0.13 

-0.94 ±0.21 ±0.12 

-0.96 ±0.42 

-1.1 ±0.4 


Table 1 

Comparison of the aA^ result with other experiments. 

4^.1 Check of the statistical errors returned by the fitter and possible bias 


A check was performed to assure that the errors returned by the fit were 
realistic and to assess if there was a bias in the fitting procedure. For both the 
particle and the antiparticle cases, the bin contents of the four cos6' region 
mass plots (Fig. 2 and 3) and the four wrong sign plots were randomly varied 
according to Poisson statistics 1000 times. Each set of eight ’’varied” plots was 
analyzed with the fitting program as it were regular data and the values of 
A and extracted. With this method we assessed a bias of 0.000 ± 0.005 
on A and we simply added 0.005 in quadrature to the systematic error. For 
oAc we assessed a bias of —0.010 ± 0.005 and a slight underestimation of 
the statistical error returned by the ht (0.16 instead of 0.15). Consequently 
this bias was subtracted, an error of 0.005 was added in quadrature to the 
systematic error and the statistical error was increased to 0.16. 


5 Final results and comparison with previous results and with the 
theoretical predictions 


After all systematic checks, bias correction and error inflation the final results 
of the measurements are = —0.78±0.16±0.13 and A = —0.07±0.19±0.12. 

Previous results obtained for are listed in Table 1. The present analysis 
result is a higher precision and it is in agreement with earlier results. No 
previous measurements for A have been reported. 

Theoretical predictions of the two-body baryon-pseudoscalar decays of charmed 
baryons have been made by a number of authors in the past. They predict the 
parity conserving amplitude A_|_ and parity violating amplitude A_ of Eqn. 
(4) from which oac ^md the decay rate can be derived. The most recent pre¬ 
dictions are listed in Table 2. All predictions are in agreement within errors 
with the FOCUS experimental result. 
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Paper 


Ivanov et al. [15] (1998) 

-0.95 

Zenczykowski [16] (1994) 

-0.99 

Zenczykowski [17] (1994) 

-0.86 

Uppal et al. [18] (1994) 

-0.85 

Cheng and Tseng [19] (1993) 

-0.95 

Cheng and Tseng [20] (1992) 

-0.96 

Xu and Kamal [21] (1992) 

-0.67 

Korner and Kramer [22] (1992) 

-0.70 


Table 2 

Summary of the most recent theoretical predictions on (under the CP conser¬ 
vation hypothesis) 

6 Conclusions 


Using data from the FOCUS (E831) experiment at Fermilab, we have studied 
the decay mode A+ —> and measured the asymmetry parameter 

(under the CP conservation hypothesis) and the CP parameter A . Our result 
for oAc is consistent with previous measurements and is more accurate. Our 
result for A is the first measurement of this quantity and is consistent with 
zero. 
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Events/ 5 MeV/c Events/ 5 MeV/c' 



Fig. 1. Avr effective mass after the selection, a) total sample; b) Avr"*" sample; c) Avr 
sample. 
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Fig. 2. Att"*" effective mass for the A+ subsample in each cos 6 bin: a) 
—1 < cos 9 < —0.5; b) —0.5 < cos 9 < 0; c) 0 < cos0 < 0.5; d) 0.5 < cos 9 < 1. 
The results of the analysis fit (solid line) are shown superimposed on the plots. The 
A^ —Att^ contribution to the fit is shown by the solid histogram, the A+ ^ 
contribution by the hatched histogram. 
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Fig. 3. Avr" effective mass for the A“ subsample in each cos 9 bin: a) 
— 1 < cos 9 < —0.5; b) —0.5 < cos9 < 0; c) 0 < cos0 < 0.5; d) 0.5 < cos9 < 1. 
The results of the analysis fit (solid line) are shown superimposed on the plots. The 
A~ ^ Avr” contribution to the fit is shown by the solid histogram, the A“ —> S 7r“ 
contribution by the hatched histogram. 




Fig. 4. Correction function e(0j, a) calculated from Monte Carlo, for the A;)' sample, 
at ten values of xcka, spanning the range from —0.9 to 0. All correction functions 
are normalized such that they are 1 in the first cos 9 bin. Note that an overall 
normalizing factor can be chosen for the correction functions without changing the 
result of the fit. 
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Fig. 5. Correction function e(0j, a) calculated from Monte Carlo, for the A“ sample, 
at ten values of a—x a^, spanning the range from —0.9 to 0. All correction functions 
are normalized such that they are 1 in the first cos 9 bin. 



cos6 


Fig. 6. Visualization of the fit result for the A+ —> A7r+ sample. The way the values 
of the solid circles and their errors are obtained are explained in the text. The solid 
straight line is 1 + q;Ac®a cos 0* where cos 6i is the value in the middle of the i**^ bin 
and OAc X Qia is the value returned by the fit. 


15 

























































acceptance corrected, normalized yield 



cos6 


Fig. 7. Visualization of the fit result for the A“ —> Att” sample. The way the values 
of solid circles and their errors are obtained are explained in the text. The solid 
straight line is 1 + a—aj^ cos 9i where cos 9i is the value in the middle of the bin 
and is the value returned by the fit. 
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Fig. 8. The correction function e(0j, a) corresponding to x oa = —0.9 (circles) 
and to a Ac xua = 0 (crosses) for the A+ sample, plotted superimposed to each other. 
In this figure a means actually cjAc x ctA- All the correction functions corresponding 
to intermediate values, for any given cos 9 interval, lie between the values plotted 
on this figure. 
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Fig. 9. The correction function e(0i, a) corresponding to a— x = —0.9 (circles) 
and to a—xaj^ = 0 (crosses) for the A“ sample, plotted superimposed to each other. 
In this figure a means actually a— x All the correction functions corresponding 
to intermediate values, for any given cos 9 interval, lie between the values plotted 
on this figure. 
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Fig. 10. Summary of the systematic checks performed on q;a^, calculated in the 
hypothesis of CP conservation. The first circle on the left of plot is the result 
quoted in this paper. 
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